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! ABSTRACT 



Aims. We carry out the effect of /3-prescription for viscosity which introduced by Duschel et al. 2000 & Hure, Richard & Zhan 
2001, in a standard self-gravitating thin disks. We were predicted in a self-gravitating thin disk the /3-model will have different 
dynamical behavior compare the well known a-prescriptions. 



■.^£2 Methods. We used self-similar methods for solving the integrated equations which govern the dynamical behavior of the thin 
f****) , disks. 

Results. We present the results of self-similar solutions of the time evolution of axisymmetric, polytropic, self-gravitating viscous 
' disks around a new born central object. We apply a /3-viscosity prescription which has been derived from rotating shear flow 
I i experiments (y = [3r 2 Q). Using reduced equations in a slow accretion limit, we demonstrate inside-out self-similar solutions 
Q ' after core formation in the center. Some physical quantities for /3-disks are determined numerically. We have compared our 
results with a-disks under the same initial conditions. It has been found that the accretion rate onto the central object for 
CZ) ■ /3-disks more than a-disks at least in the outer regions where /3-disks are more efficient. Our results show that Toomre instability 
parameter is less than one everywhere on the /3-disk which means that in such disks gravitational instabilities can be occurred, 
so the /3-disk model can be a good candidate for the origin of planetary systems. Our results show that the /3-disks will de- 
couple in the outer part of the disk where the self-gravity plays an important role which is in agreement with Duschl predictions. 
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1. Introduction 

Accretion disks are recognized as one of the objects that 
are found around many astrophysical objects, such as ac- 
tive galactic nuclei AGN, binary stars and young stellar 
objects. On the observational side, evidence for disks in 
young stellar objects gleaned both spectroscopically and 
through direct imaging is now quite compelling (Beckwith 
& Sargent 1993, Storm et al, 1993). Up to half of the solar 
type, pre-main sequence stars are surrounding by diks of 
gas and dust, many of these disks having masses similar to 
that expected for the early solar nebula (Chandler, C. J. 
1998). The accretion disks around pre-main sequence stars 
are good candidates for the creation of planetary systems. 
The structure of such disks is a subject of great interest 
and has been studied both through self-similar solutions 
assuming unsteady state (Mineshige & Umemura 1997; 
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Mineshige et al. 1997; Tsuribe 1999) and through direct 
numerical hydrodynamical simulations (Igumenshchev & 
Abramowicz 1999; Stone et al. 1999; Torkelsson et al. 
2000). It is understood that the most crucial factors are 
self-gravity and viscosity which have great role on angu- 
lar momentum transportation on the gas disk. Accretion 
takes place because of the action of some form of dissi- 
pation which releases the free energy of the shear flow as 
heat, and so allows the disk material to fall deeper into 
the potential well of the central object. In a simple pic- 
ture ,Lynden-Bell and Pringle (1974) indicated that the 
dissipative processes must take the form of a stress which 
transport angular momentum outwards. It plays a signifi- 
cant role in many such systems, ranging from protostellar 
disks to active galactic nuclei (AGN). The self-gravity will 
modify the radial and vertical equations and so can treat 
the dynamical behavior of the accretion disks. In the stan- 
dard thin accretion disk model, the effect of self-gravity is 
neglected, and only pressure supports the vertical struc- 
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turc. By contrast, the theory of self-gravitating accretion 
disks is less developed and in the traditional model of ac- 
cretion disks, the self-gravity is ignored just for simplicity 
(e.g. Pringle 1981), although the self-gravity can describe 
the deviation from Keplerian rotation velocity in some 
AGN and flat infrared spectrum of some T Tauri stars. 
From the observational point of view, there are already 
some clues that the disk self-gravity can be important 
both in the context of proto-stellar disks and in the accre- 
tion disks around super massive black holes in the AGN. 
However, the detailed comparison with observations is lim- 
ited by the lack of detailed models of self-gravitating disks 
and by an incomplete understanding of the basic physical 
processes involved. 

The study of self-gravity in a most general case is dif- 
ficult and because of these complexities, authors usually 
study the effects related to the self-gravity either in the 
vertical structure of the disk (e.g. Bardou et al. 1998) or in 
the radial direction (e.g. Bodo & Curir 1992). Disks in the 
AGN are thought to be relatively light in the sense that 
the ratio of is around a few percents (where M and M* 
are the masses of the the disk and central star). Usually 
self-gravity plays its role at large distances from the cen- 
tral objects (Shlosman & Begelman 1987), and mainly in 
the direction perpendicular to the plane of the disk. But 
in the accretion disks around young stellar objects or pre- 
main-sequence stars, self-gravity can be important in all 
parts of the disk in both vertical and radial directions. 
Early numerical works of self-gravitating accretion disks 
began with N-body modelling (Cassen & Moosman 1981 
; Tomley, Cassen & Steinman- Cameron 1991). Shlosman 
& Begelman (1987) investigated the role of self-gravity in 
AGN. Recently, Ghanbari and Abbassi (2004) introduced 
a toy model which shows that at least in equilibrium of a 
thick accretion disk, self-gravity is an important effect. 

One of the basic concepts of the theoretical descrip- 
tions of accretion disks and these dynamics, is the de- 
tailed knowledge underlying the physics of viscosity in the 
disks. Because of all detailed modeling of the structure 
and evolution of accretion disks depending on the viscos- 
ity and its dependence on the physical parameters, choos- 
ing the best viscosity model is quite important. There is 
a belief that, the molecular viscosity is inadequate to de- 
scribe the observational evidence of some luminous accre- 
tion disks so some kinds of turbulence viscosity are re- 
quired. Most investigators adopt so-called a-model intro- 
duced by Shakura (1972) and Shakura & Sanyeav (1973) 
that gives the viscosity as the product of pressure scale 
height in the disk (h), the velocity of the sound (c s ) and 
a parameter a which contains all unknown physics. The 
models for the structure and evolution of accretion disks 
in close binary systems (e.g. dwarf novae and symbiotic 
stars) show that Shakura and Sunyeav's parametrization 
with a constant a leads to results that reproduce the over- 
all observed behavior of the disks quite well. And it has 
been recently recognized that accretion disks treated by 
a week magnetic field arc subjected to MHD instabilities 
(Balbus & Hawly 1991), that can induce some kind of 



turbulence in the disk, thereby being able to transport 
angular momentum and to promote accretion processes. 
However in many astrophysically interesting cases, such 
as the outer part of proto-stellar disks, the ionization level 
is expected to be low, reducing significantly the effect of 
magnetic field in the dynamical behaviour of the disk. The 
realization that molecular transport of angular momen- 
tum is so inefficient led the theoreticians to look for an- 
other mechanism of transport of angular momentum in 
accretion disks. The alternative mechanism that would be 
a good candidate for transporting angular momentum in 
accretion disks is any kind of turbulence. Actually the 
a-prescription is based on a kind of turbulence viscosity 
but it has not any physical base to drive the origin of 
turbulence in the model. On the other hand, laboratory 
experiments of Taylor-Couette systems seem to indicate 
that, although Coriolis force delays the onset of turbu- 
lence, the flow is ultimately unstable to turbulence for 
Reynolds numbers larger than a few thousands (e.g., see , 
Richard, Zahn, 1999 and Hure, Richard, Zahn 2001). Since 
in all kind of self-gravitating disks the Reynolds number 
is extremely high, it was thought that the hydrodynam- 
ical driven turbulent viscosity based on critical Reynolds 
number has probably significant role in the distribution 
of angular momentum in the accretion disks. The result- 
ing turbulence would then transport angular momentum 
efficiently. Recently, Duschl, Stirittmatter and Bicrmann 
(2000) have proposed a generalized accretion disk viscosity 
prescription based on the hydrodynamically driven turbu- 
lence at the critical effective Reynolds number,/?- mo del , 
which is applied for both self-gravitating and non self- 
gravitating disks and is shown to yield the standard a- 
model in the case of the shock dissipation limited, non 
sclf-gravitating disks. They have shown that in the case 
of fully self-gravitating disks this model may explain the 
observed spectra of proto-planetary disks and yield a nat- 
ural explanation for the radial motions from the observed 
metallicity gradients in the disk galaxy. 

The basic equilibrium and dynamical structure of ac- 
cretion disks are now well understood, as long as the stan- 
dard model based on the a-viscosity prescription (Shakura 
& Sunyeav 1973) is believed. Nevertheless, it is not easy 
to follow its dynamical evolution, mainly because the ba- 
sic equations of the system are highly non-linear, spe- 
cially when the system is self-gravitating (Paczynski 1978; 
Fukuc & Sakamoto 1992). To follow the non-linear evolu- 
tion of dynamically evolving systems, in general, the tech- 
nique of the self-similarity is sometimes useful. Self-similar 
assumptions enable us to simplify the governing equations. 
Self-similar solutions have a wide range of applications in 
astrophysics. Several classes of self-similar solutions were 
known previously (Pringle 1974, Filipov 1984), but all of 
them considered the disk in a fixed, central potential. But 
a class of self-similar solution had been provided dur- 
ing last years that contained self-gravity (Mineshige & 
Umemura 1997). They had found a self-similar solution for 
a time evolution of isothermal, self-gravitating a-viscous 
disk. This solution describes a homologous collapse of a 
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disk via self-gravity and viscosity. They found that the 
disk structure and evolution are distinct in the inner and 
outer parts. The effect of self-gravity in the collapse of 
polytropic self-gravitating viscous disk has been investi- 
gated by Mineshige, Nakayama and Umemura (1997). 

Following the Duschl et al. (2000) suggestion for (3- 
prescription for viscosity, we apply this model for a thin 
self-gravitating disk around newborn stars. At first, it may 
seem that using other forms of viscosity is not an im- 
portant issue, because just one should change the math- 
ematical forms of the equations. But this simplified ap- 
pearance did not force authors from such studies because 
of possible interesting results that one may obtain from 
other viscosity prescriptions. All these lead us to explore 
a self-gravitating disk using other viscosity prescriptions. 
However, all viscosity prescriptions are suffering from hav- 
ing phenomenological backgrounds rather than physically 
confirmed backgrounds. We think while we don't have a 
clear picture of the turbulence in disks, all such prescrip- 
tions are standing on the same level as for their physical 
backgrounds. On the other hand, self-gravitation in a disk 
is a highly nonlinear process as a result of the complex 
behaviors of the various physical agents of the system, in 
which the turbulent viscosity and its prescriptions have a 
vital role. Thus, one may naturally ask what would happen 
if other forms of viscosity prescriptions in a self-gravitating 
disk is being used. When we searched in the literature, 
we found the /3-prescription as only experimentally tested 
viscosity prescription. Indeed, there are few studies about 
/3-prescription comparing the a-model, clearly because of 
its newly proposed form. One may refer to many papers 
to see physical considerations and applications of the (3- 
model: Mayer & Duschl (2005); Weigelt et al. (2004); Pott 
et al. (2004); Granato et al. (2004); Mathis, Palacios & 
Zahn (2004); Richard & Davis (2004). We expect to find 
different dynamical behaviors and we will show that in 
these disks the gravitational fragmentation can take place 
everywhere in the disks. So it will be a good description 
for the formation of a proto-planetary disk. 

2. Formulation of Equations for Self-Similar 
Variables 

Self-similar behavior provides a set of unsteady solutions 
to the self-gravitating fluid equations. On the other hand, 
many physical problems often attain self-similar limits for 
a wide range of initial conditions. Also self-similar prop- 
erties allow us to investigate properties of the solutions in 
arbitrary detail, without any of the associated difficulties 
of numerical hydrodynamics. 



disks are geometrically thin in the vertical direction 
and symmetric in the azimuthal direction. The model is 
described by the fundamental governing equations which 
are written as follow : 



d P , 1 9 I 2 x n 



dv r 



dv r 
dr 



1 dp 
p dr 



d<f> 
dr 



pv r 



dr 



1 



_d_ ( 4fM) 
r 2 dr dr 



(1) 



(2) 



(3) 



where p, p, v r and v$ are density, pressure, radial and 
azimuthal components of velocity of the gas disk and also 
$ is the gravitational potential of the gas disk inside of 
the radius r. We assume a polytropic relation between the 
gas pressure and density : 



P = Kp 1 



(4) 



with K and 7 being constants. The polytropic index 7 
describes the adiabatic pressure-density relation. In sub- 
sequent analysis, we vary it and represent its effect on 
some physical quantities. The vertical extent of the disk 
at any radius is given by h, the half thickness of the disk: 



h = 



(4irGp)2 



2ttGct 



(5) 



where c s is the sound speed and a is the surface density. 

The solution of these equations, give us the dynam- 
ical evolution of the disk which strongly depend on the 
viscosity model. So the study of dynamical behavior of 
the accretion disks is postponed to the more information 
about the viscosity. 

2.2. Nondimensionalization 

Before we actually begin solving equations Q-©, it 
is convenient to nondimensionalize the equations. The 
essence of self-similar model is the existence of only two 
dimensional parameters in the problem, viz., K and 
the gravitational constant G. It is assumed, and this is 
born out by numerical calculations, that any additional 
parameters, such as the initial central density affects only 
transients, and theis memory is quickly lost, at least in 
the part of the flow in which the density greatly exceeds 
the initial central density. If that is the case, then only 
one dimensionless combination of radius r and time t can 
be found (Mineshige et al. 1997, Yahil 1983): 



2.1. The Basic Equations 

In order to study the accretion processes of a thin disk 
under the effect of the self-gravity and viscosity, we 
consider axisymmetric polytropic disks using the cylin- 
drical coordinates (r,(j),z). We assume that the accretion 



i = K- 1 2G Z ^r{ty- 2 (6) 

This determines the dimensionless parameterizations 
of any similarity solution (Mineshige et al. 1997). Note 
that, we consider only t > for our work and the origin 
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of time (t = 0) corresponds to the core formation epoch. 
Hence we have: 



d_ 

dt> 



(7- 



(7) 



d_ 

at ~ 

and 

l=K-^( t r-^ (8) 

for the transformation (t, r) — > (t' — t,£). Self-similarity 
allows us to reduce the self-gravitating fluid equations 
from partial differential equations into ordinary differen- 
tial equations. 

For changing the variables to dimensionless form, we 
used K and G; because we require that all of the time- 
dependent terms disappear in the self-similar forms of the 
equations. Other physical quantities (functions of t and r) 
are transformed into self-similar ones (functions of only £) 
as: 



v T {t,r) = K^G^{t) x -^V r {i), 
v 4> (t,r)=KiG^ L (t) 1 -"V lP (a 

a(t,r) = (27T)- 1 K^G-^ L (t)-^), 
p(t,r) = (47T 7 )-* K*G-\t)- 2 tf (0, 
p(t,r) = (An^KG-^ty^iO, 
KG x -^{tf-^v'{0 

KiG^itf-^M^). 



u{t,r) = 
m(t, r) = 
where 



M 



(9) 
(10) 
(11) 
(12) 
(13) 
(14) 
(15) 
(16) 

(17) 



3 7 -4 

with respect to another form for the continuity equation 



dm dm 
dt r dr 

thus, we have 







m = 2tt 



ardr 



M 



3 7 -4 



(18) 



(19) 



In order to solve the equations, we need to assign the 
kinematic coefficient of viscosity v. Although there are 
many uncertainties about the exact form of viscosity, as we 
mentioned in Introduction, we employ the /3-prescription 
introduced by Duschl et al.: 



(20) 



where v' is on a dimensionless form. We expect when 
the disk is non-self-gravitating, it leads to standard a- 
prescriptions. With substituting this prescription into the 
above equations we can investigate the dynamical evolu- 
tion of the disk. 



2.3. Transformation of the basic Equations 

Substituting the above transformations in equations fT|l- 
© , we can introduce a set of coupled ordinary differential 
equations. The basic equations are then transformed into 
the following forms: 



1 d 



(£E«) = (3 7 - 4)S 



(21) 



du c 2 dS M J 2 s w ^, , 



d J 1 d , o d , J „ 
where u is defined as, 



u = V r -(2- 7 )£ 
and c 2 is a constant: 



1 — 7 27-2 

c 2 = 2(Attj)~T,^T 



(23) 



(24) 



(25) 



Now we have a set of complicated differential equa- 
tions that must be solved under appropriate boundary 
conditions. Although a full numerical solutions to these 
equations would now be possible; it is more instructive to 
proceed by analyzing the model in some restrictive cases 
such as one on slow accretion limits. 

3. Behavior of the Solutions 
3.1. Slow Accretion Limit 

We consider the fluid equations for a thin disk in the slow 
accretion approximation. The most substantive aspect of 
the slow accretion approximation consists of rotationally 
supported disks when the viscous timescale is much larger 
than the dynamical timescale. In addition, the pressure 
gradient force and the acceleration term in this approxi- 
mation are ignored. The slow accretion approximation in 
disks has been used by Tsuribe (1999) and Minshige et al. 
(1997) and many others. 

In the slow accretion limit (\V r \ <C 1 -C V^), and in 
the equation |3 the Euler equation, we have radial force- 
balance which means that only two terms on the right 
hand of the equation balance each other so we have: 







J 2 M 
hence we have: 



J=(£M)i = ( 



3 7 -4 ; 



(26) 



(27) 



We can make some simplifications in the equation of 
continuity in self-similar form, then it becomes: 



dlnY] 
dln£ 



-1 



dlnu 37 
dlnl + (_ i 



(28) 
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Before driving other equations in a slow accretion 
limit, we introduce the distribution of the initial specific 
angular momentum, j, as (Tsuribe 1999): 



G 

j(r) — q — m(r) 



(29) 



where j and m are self-similar angular momentum and 
the total mass within the cylindrical radius r, respectively. 
The variable q is a dimcnsionless quantity has a constant 
value for the invicid equilibrium solutions (Mestel 1963; 
Toomrc 1982) and variable value when the viscosity effect 
is included. Thus in dimcnsionless form, 



J = qM 

Now using Eas. (j2*7jl and (|3T)|> we find that: 
(37-4)4 



qu< 



(30) 



(31) 



With Eqs. (ESI,®,® .l|2g|l and (|5Tj) we can finally 
find a simple differential equation for V r : 

dV = g K[KT(2-7K] 2 [3(2- 7 )±5(3 7 -4)]K 
d£ (3 S[±3V r T (37 - 2)f] [±3V r T (37 - 2)£] 



[(2- 7 )(-18 7 + 22)T2(3 7 -4) 2 ]e 
[±W r T (37 - 2)£] 



(32) 



Now we can use the standard fourth-order Runge- 
Kutta scheme to integrate this ordinary differential equa- 
tion. We can also use the asymptotic solutions of this 
equation near the origin and the outer part of the disk 
as a boundary condition. 

3.2. Numerical Analysis 

To solve equation l |3"2"|l numerically, we need one boundary 
condition. We derived asymptotic solutions in the limits 
£ — ► and £ — > oo for V r where these asymptotic values 
can be used as boundary condition. 



V r = 



5 7 — 6 



<?o (7 - I) 2 



6 7 -7 L " 13 6 7 -7 
[(2 - 7 )(22- 18 7 )- 2(3 7 -4) 2 ] 



(2-7) 5 



(33) 



£ ^ oo (34) 



where go is the asymptotic value of q in the outer ra- 
dius. The equation (|32() is integrated in the limits of the 
following equations (eqs. 33 and 34). It was found that 
only for some special values of the physical parameter 
such as <7, 7 ,/3 there exist physical solutions satisfied at 
both boundary conditions. By solving this equation we 
will have V r profile as a function of £ and other profiles 
(V0,£M...) can be obtained easily. In figures 1 and 2, we 
show radial and azimuthal velocity distributions for some 
7 and (3 values, respectively. Figures 3 and 4 indicate sur- 
face densities and mass accretion rates profiles for some 7 
and (3 values. Also we compare (3 and a disks for 7 = 1.1 



in all figures. The self-similar variables are functions of 
£. The behavior of the solutions predicted by the beta 
viscosity makes much more radial velocity (see Figure 1) 
compared to the a-model (at least in the outer part of 
the disk where the self-gravity has an important role) and 
also for the azimuthal component of self-similar velocity, 
/3-viscosity leads to faster velocities compared to the a- 
model (see Figure 2). Therefore, the viscosity has more 
efficient role on the redistribution of angular momentum, 
and it leads to more radial flow and accretion rates. 

It is predicted that the outer parts of a thin, viscous 
disk around QSOs, self-gravity has an important role. This 
effect is investigated by Toomre parameter (Toomre 1964), 
such that the local gravitational instability occurs where 
Q < 1 and where Q > 1 the disk is stable against the 
gravitational fragmentation. So it is useful to calculate Q 
values to compare gravitational stability in (3 and a disks. 
The Toomre stability parameter for an epicyclic motion 



Q = -pr = L, y K— 4tt — 7 2^s~ 
7tG(7 tr 

where 

k = 0(4 + 2^)4 
dlogr 



(35) 



(36) 



is the epicyclic frequency and j£ — jy. So if 7 = 1, we 

obtain Q = (Tsuribe 1999). 

In figure 5, we show the distribution of the Toomre Q 
value for some 7 and (3 values. We also compare (3 and 
a disks for 7 = 1.1. In figure 6, the angular momentum 
coefficient q = j-f is plotted as a function of £ for some 7 
values at j3 = 10 -3 . Also we see its behavior for (3 and a 
disks. 



4. Time scales 

For estimating the effect of viscosity on the evolution of 
accretion disks we can compare the viscose time scale with 
dynamical time scale. The dynamical time scale Td yn is 
given by: 



1 

Tdyn - 



(37) 



where f2 in the non self-gravitating (NSG) and Keplerian 
self-gravitating (KSG) is given by the mass of the central 
accretor and by the radius. But in the fully self-gravitating 
disks, O is determined by solving Poisson's equation. The 
time scale of viscose evolution r V i SC is given by: 



r 
v 



(38) 



where in the standard non self-gravitating and geometri- 
cally thin accretion disks where h « r ( case a-disks ), 
this leads to 



non — SG 



1 r \ 2 T dyr, 

[ h> a 



(39) 



6 



Shahram Abbassi et al.: Self-Similar Evolutionary Solutions of Self-Gravitating, Polytropic j3- Viscous Disks 



where in full self-gravitating (FSG) or KSG disks ( (3- 
disks) is given by 



FSG _ SG _ KSG 
vise vise vise 



T~dyn 



(40) 



So with a < 1 and p << 1 under all circumstances 
T v isc >> i~dy n - So the slow accretion limit will be con- 
firmed. The best approximation for a parameter in the 
standard accretion model is less than one, a ~ 0.1 but in 
the [3 model which based on the critical Reynolds number 
beta is approximately (3 ~ 10~ 2 — 10 -3 . (-) in the outer 
part of accretion disks is too small, so in the outer part, of 
the disk where self-gravity of the disk dominates the beta 
viscous time scale is less than alpha viscous time scale. 
Therefore, the dynamical evolution of /3-disks at least in 
the outer part is faster than the inner parts. Although, 
in the case of non-SG disks /3-model can recover the stan- 
dard a-disks. So with /3-model we can reconstruct a better 
picture in the equilibrium of galactic disks and protoplan- 
etary disks. As we can see in the Fig.l, in comparison to 
alpha model, the radial flow in the beta disk is quite high 
which relevant to beta mechanism time scale. Because of 
the time scale of the beta disk is less, so it can evolve 
disks effectively and it can produce the large radial flow. 
In the case of galactic disks the inflow velocities in the 
/3-model suggests values in the range 0.3 — 3fcms _1 which 
is quite difficult to measure directly where a-model sug- 
gest still lower values Duschl et al. (2000). Many authors 
have suggested the radial abundance gradients observed in 
our own and other disk galaxies maybe due to radial mo- 
tion and diffusive mixing associated with the turbulence 
generated by eddy viscosity ( /3-disks) (Lacey & Fall 1985, 
Sommer-Larsen & Yoshii 1990; Koppen 1994; Edmunds & 
Greenhow 1995). Such radial inflows are consistant with 
the /3-model, but could be described with other physical 
processes such as the effect of magnetic field. 

5. Concluding Remarks 

The /3-prescription based on the assumption that the ef- 
fective Reynolds number of the turbulence dose not fall 
bellow the critical Reynolds number. In this parametriza- 
tion the viscosity is proportional to the azimuthal veloc- 
ity and the radius. This model yields physically consis- 
tent models of both Keplerian and fully self-gravitating 
accretion disks where in the case of thin disks with suffi- 
ciently small mass, recover the a-disk solutions. Such (3- 
disk models may be relevant to protoplanetary accretion 
disks (Duschl et al. (2000)). In the case of protoplane- 
tary disks they yield spectra which are considerably flat- 
ter than those due to non-self-gravitating disks, in better 
agreement with observed spectra of these objects. 

In this paper, we have considered the time-dependent 
evolution of self-gravitating disks with /3-prescription by 
self-similar method for a thin, viscous disk. To do this, 
we started from dimensionless basic fluid equations. In 
order to dominate gravity and the centrifugal force, we 
consider the fluid equations for a thin disk in the slow 




Fig. 1. Radial self-similar velocity distributions as a func- 
tion of self-similar variable £ for a) 7 = 1.0,1.1 at /3 = 
10~ 3 b) p = 10- 3 ,2xl0- 3 ,3xl0~ 3 at 7 = 1.1 c) (3 = 10~ 3 
and a = 10 _1 at 7 = 1.1. 



2 -1 
-2 
1 

„" 
> 

S -1 



> -2 

§> 
o 

-4 















P=1(T 3 












H-1 


(a) 


3 


-2 


-1 


1 2 

log® 


3 


4 


5 6 














-,.1.1 














lb) 


3 


-2 


-1 


1 2 

log(5) 


3 


4 


5 6 












f> 


-,.1.1 




























(c) 



-3 -2 -1 1 2 3 4 5 6 

logs 

Fig. 2. Azimuthal self-similar velocity distributions as a 
function of self-similar variable £, corresponding to a) 7 = 
1.0,1.1 at (3 = 10- 3 b) (3 = 10~ 3 ,2 x 10~ 3 ,3 x 10~ 3 at 
7 = 1.1 c) $ = 10~ 3 and a = 10" 1 at 7 = 1.1. 



accretion approximation. It has been found that evolution 
is described by solving a simple differential equation (|32() . 
We solved it numerically, beginning asymptotic solution 
of this equation near the origin as a boundary condition. 
Note should be taken that, we had the limitation to select 
parameter 7 for essence of differential equations and the 
fact that we seek just physical solutions. 

The presented method shows that an increase of P 
value causes the azimuthal velocity to decrease but its 
general distribution function doesn't vary throughout the 
disk. Also, azimuthal velocities in /3-disks are far more 
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Fig. 4. Self-similar mass accretion rates distributions as 
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Fig. 5. Distribution of the Toomrc Q value as a function 
of self-similar variable £, corresponding to a) 7 = 1.0, 1.1 
at /3 = 10- 3 b) 13 = 10~ 3 , 2 x 10~ 3 , 3 x 10~ 3 at 7 = 1.1 c) 
13 = 10~ 3 and a = 10" 1 at 7 = 1.1. 
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Fig. 6. The angular momentum coefficient profile, q, as 
a function self-similar variable £ corresponding to a) 7 = 
1.0,1.1 at (3 = 10~ 3 . We see its behavior versus the di- 
mensionless viscosity parameter for b) (3 disks c) a disks 
at 7 = 1.1. 



than a-disks (see Figure 2). So we expect the /3-disks 
evolve in different ways with respect to the a-disks. 

According to Figure 4, /3-disks have larger mass accre- 
tion rates than a-disks. So, observably, we expect them 
to be brighter than a-disks. Also, we note that with the 
increase of the (3 value, M increases. Then mass flows in- 
crease onto the central object. This follows more radial 
velocity and less surface density (see Figures 1,3). 

Comparing to a-disks, q distribution (q = -jg) seems 
non-smoothly. The q values are very small in the inner- 
most regions (see Figure 6). Whereas it is almost constant 
in the outer regions. It seems in the outer part of the 



disk where the beta viscosity is more efficient, the angu- 
lar momentum is proportional to the disk's mass inside 
the radius r. In order to study the effect of self-gravity 
of thin /3-disks, we plot the Toomre parameter as a func- 
tion of £. It is obvious that the gravitational instabilities 
in /3-disks are more pronounced than a-disks, Q < 1. In 
Figure 5, Toomre parameter profile can reveal this sub- 
ject. So it can be expected that the /3-disk is a good model 
to describe planet formation around new-born stars. In a 
global overview, we showed that in the outer parts of the 
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disk there is a difference between a and (3 models. These 
results were predicated by Duschl et. al (2000). 

Further, in order to study actual model and make a 
realistic picture for a thin self-gravitating disks, one must 
investigate energy exchange of the disk with its environ- 
ment. In this case, one should find a mechanism for trans- 
ferring the thermal energy from disk to outside environ- 
ment; so we should add energy equation to our model. 
Both a and [3 models are phenomenological prescriptions 
for disk viscosity. In an actual model of viscosity, it is pos- 
sible to combine these two models and establish an exact 
description for different regions of the disk. Also in real 
accretion disks, there are many important processes other 
than viscosity that we have neglected. For example, non- 
axisymmetric waves which are also expected to transport 
angular momentum outward. Magnetic field and its in- 
fluence arc neglected and sometimes magnetic braking is 
another possibility for transporting angular momentum. 
However, these preliminary solutions can be the begin- 
ning of our understanding deployment from the physics 
governing the accretion disk. 
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